Stiffness in stochastic chemically reacting systems: The implicit tau-leaping method
نویسندگان
چکیده
We show how stiffness manifests itself in the simulation of chemical reactions at both the continuous-deterministic level and the discrete-stochastic level. Existing discrete stochastic simulation methods, such as the stochastic simulation algorithm and the ~explicit! tau-leaping method, are both exceedingly slow for such systems. We propose an implicit tau-leaping method that can take much larger time steps for many of these problems. © 2003 American Institute of Physics. @DOI: 10.1063/1.1627296#
منابع مشابه
The numerical stability of leaping methods for stochastic simulation of chemically reacting systems.
Tau-leaping methods have recently been proposed for the acceleration of discrete stochastic simulation of chemically reacting systems. This paper considers the numerical stability of these methods. The concept of stochastic absolute stability is defined, discussed, and applied to the following leaping methods: the explicit tau, implicit tau, and trapezoidal tau.
متن کاملExposition and Streamlined Formulation of Adaptive Explicit-Implicit Tau-Leaping∗
The adaptive explicit-implicit tau-leaping method with automatic tau selection [Y. Cao, D. T. Gillespie and L. R. Petzold, J. Chem. Phys., 126(22):224101 (2007)] is a flexible algorithm for accelerated stochastic simulation of chemically reacting systems. It combines the advantages of different simulation schemes and is particularly useful when a system changes its dynamical behavior over time ...
متن کاملTau leaping of stiff stochastic chemical systems via local central limit approximation
Stiffness manifests in stochastic dynamic systems in a more complex manner than in deterministic systems; it is not only important for a time-steppingmethod to remain stable but it is also important for the method to capture the asymptotic variances accurately. In the context of stochastic chemical systems, time stepping methods are known as tau leaping. Well known existing tau leaping methods ...
متن کاملComment on "Stiffness in stochastic chemically reacting systems: the implicit tau-leaping method" [J. Chem. Phys. 119, 12784 (2003)].
sions can be obtained by using the generating function in solving the master equation. In addition, our approach allows for the computation of the probabilities for the number of S1 molecules at a certain time t when the reaction parameters, i.e., the initial number of molecules X1(0) and X2(0) of S1 and S2 , respectively, and the reaction rates c1 and c2 , are given. The master equation is a d...
متن کاملThe Adaptive Explicit-Implicit Tau-Leaping Method with Automatic Tau Selection
The existing tau-selection strategy, which was designed for explicit tau-leaping, is here modified to apply to implicit tau-leaping, allowing for longer steps when the system is stiff. Further, an adaptive strategy is proposed that identifies stiffness and automatically chooses between the explicit and the (new) implicit tau-selection methods to achieve better efficiency. Numerical testing demo...
متن کامل